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We study the azimuthal angular dependence of back-to-back di-jet production in unpolarized 
hadron scattering Ha + Hb — > Ji + Ji + X, arising from the product of two Boer- Mulders functions, 
which describe the transverse spin distribution of quarks inside an unpolarized hadron. We find that 
when the di-jet is of two identical quarks (J q + J q ) or a quark-antiquark pair (J q + J q ), there is a 
cos 8(f) angular dependence of the di-jet, with S<j> = 4>i — 4>2 , and (f>i and (j>2 are the azimuthal angles 
of the two individual jets. In the case of J q + J q production, we find that there is a color factor 
enhancement in the gluonic cross-section, compared with the result from the standard generalized 
parton model. We estimate the cos<5</> asymmetry of di-jet production at RHIC, showing that the 
color factor enhancement in the angular dependent of J q + J q production will reverse the sign of the 
asymmetry. 

PACS numbers: 12.38.Bx, 13.85.Ni, 13.88.+C 

I. INTRODUCTION 

The role of transverse momentum dependent (TMD) distributions [H, i n high energy physics has received a great 
deal of attection recently. This major interest in TMD distributions relies on the belief that they are responsible for 
various azimuthal asymmetries observed in different single particle inclusive processes involving at least two hadrons. 
Examples are the single transverse spin asymmetry (SSA) measured in semi-inclusive deeply inelastic scattering 
(SIDIS) 0,01 and inclusive pionproduction in hadron collision [1, @], as well as large cos2(/> anomalous asymmetry 
in the Drell-Yan process 0, H, It is found that leading- twist time-reversal-odd (T-odd) distributions play an 
essential role in these asymmetries, i.e. the Sivers function [l(| can account for the SSA, while the Boer-Mulders 
function 0, HH can explain the large cos 2(j> anomaly. A key ingredient of T-odd distributions is the path-ordered 
exponential contained in their definition, the so called Wilson line or gauge-link [l2|, [HI [HI ^ ue to initial-state or 
final-state interactions [To] ] between the active partons and the spectator system. For example, final state interactions 
are related to a future pointing Wilson line (denoted as U^) for T-odd distributions, connected to the SSA in 
SIDIS processes. On the other hand, the study of Drell-Yan processes shows that initial-state interactions give rise 
to a past-pointing Wilson line (U^), which predicts a minus sign difference of the Sivers functions (l3l. [lil Tl5j and 
correspondingly the SSA in the Drell-Yan process compared to the SIDIS process. Based on these studies, TMD 
Factorization in SIDIS, Drell-Yan and e + e~ annihilation has been established , [TtJ , within the allowance of a sign 
reversal for T-odd functions in different processes. 

Recently further studies [ill QjJ [2(| [2l|, H2, [H, [24|, [2|| have been performed on back to back di-jet or di-hadron 
production in hadron-hadron collisions. Unlike in SIDIS and Drell-Yan processes, in which there is only either final- 
state interactions or initial state interactions, in hadron-hadron collisions both initial-state and final-state interactions 
contribute to the process. The result is that there is a more complicated form of the gauge-link for TMD distribution 
functions in hadron-hadron collisions than in SIDIS or Drell-Yan processes. For example the loop structure = 
WH'wW will appear in the corresponding gauge- link [llj] . Further studies [2(3] show that in hadron-hadron collisions 
the form of the gauge-link depends even on the hard partonic process. The impact of these nontrivial gauge-links on 
phenomenology, especially on SSA, has been studied in back-to-back hadron pair or di-jet production [l3. l2ll. |2(| . and 
photon-jet productions [2lJ- In Refs. [H, the authors show that in second order of gluon exchange from initial- 
and final-state interaction, the standard TMD factorization is violated in back-to-back hadron (jet) pair production, 
i.e., the TMD distributions need to be modified in a process-dependent way (See also Ref. [29| , where the universal- 
breaking part of the TMD distributions has been given). 

As shown in Refs. [H, a consequence of nontrivial gauge-links is that in hadron-hadron scattering the parton 
distribution (and fragmentation) functions are convoluted with gluonic pole cross-sections rather than the standard 
basic parton cross-sections. Compared to the basic parton cross-sections, there is a gluonic pole factor Cq along 
with the gluonic cross-section, which arises from the T-odd effect contributed by the multiple initial- and final-state 
interactions mediated by gluon exchanges. The complete gluonic pole cross-section encountered in SSA in hadron- 
hadron scattering have been given in Ref. [Ill [HI] for hadron pair, di-jet and hadron-jet production, and in Ref. [13] for 
photon-jet production. In this paper, we will study the azimuthal angle dependence of back-to-back di-jet produced 
in unpolarized hadron-hadron scattering. Like the SSA in di-jet production, where one incident hadron is transversely 
polarized, there is also interesting phenomenology in the case of di-jet production in unpolarized hadron-hadron 
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FIG. 1: partonic diagrams contributing to di-jet production Ha + Hb — > Jq + Jq + X. 



scattering. First, one realizes that T-odd TMD distributions, especially the Boer-Mulders function, produces also 
azimuthal asymmetries in unpolarized hadron scattering, similar to the case of the cos 2<fi asymmetry in unpolarizcd 
Drcll-Yan processes. A recent example is the analysis of that asymmetry of photon-jet production in unpolarizcd 
hadron scattering, given in Ref. (30| . Therefore the study will provide more knowledge on the TMD distributions 
which play a role in the process. Second, it is interesting to investigate the effect of multiple initial- and final- 
state interactions on the cross-section of di-jet production in unpolarized hadron scattering, as has been done in single 
transversely polarized hadron scattering. In unpolarized hadron scattering the azimuthal asymmetry can be produced 
by two T-odd distributions, which will make the gluonic pole factor more complicated. Therefore a quantitative study 
of the effect of multiple initial- and final-state interaction on unpolarized hadron scattering will shed light on the 
QCD dynamics that is present in azimuthal asymmetries. 



II. ANGULAR DEPENDENCE OF DI-JET PRODUCTION IN UNPOLARIZED HADRON 

SCATTERING 



The process we study is 

H A (P A ) + H B (P B ) - J X {P X ) + J 2 (P 2 ) + X, (1) 

where the two incident hadrons Ha and Hb are unpolarized. The momenta of these initial hadrons are denoted by 
Pa and Pa, and those of the di-jet by Pi and P%. We are interested in the kinematical region where the transverse 
momenta of the di-jet, Pu_ and P2±, have similar size and are almost back to back 

PlX~-P2±,\Pl±\~\P2±\~P±, (2) 

and where the total transverse momentum of the di-jet q± = Pu_ + P 2 ± is much smaller than P±: q± = | q±\ -c Pi- 
Therefore the process is sensitive to the intrinsic transverse momenta of the partons inside the hadron. We label the 
azimuthal angles of P±± and P 2 _l as <j>\ and cj) 2 , respectively. We also define 5<p — tt ~ (<f> 2 — </>i), which describes the 
imbalance of the di-jet from a back to back configuration. 

The cross-section of di-jet production can be calculated from the general form 



±|M| 2 ^ 

2s 1 1 (2irf2E Pl {2irf2E P2 



d^=-\M\ \ n ' - 2 , (3) 



where the amplitude square is expressed as a convolution of the parton-parton to di-jet hard amplitudes and the 
correlation functions 11911: 



\M\ 2 = J dx 1 d 2 k 1± dx 2 d 2 k 2 ±(2TT) 4 



x <5 4 (fc! +k 2 -P 1 - P2)Tr{<S>{x 1 ,k 2 1± )$(x2,kl 1 _) 

X H(k 1 ,k2,Pl,P2)H*{k 1 ,k2,Pl,P 2 )}. (4) 

In the above equation we denote the momenta of the initial partons as k\ and fc 2 , respectively. 
The Mandelstam variables of the partonic process are defined as 

a = (h+k 2 ) 2 , (5) 
t = (fci - Pi) 2 , (6) 
u = (fci -P 2 ) 2 , (7) 
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FIG. 2: Partonic diagrams contributing to di-jet production Ha + Hb — » Jq + J q + X through qq — > qq. 



which satisfy the relation 
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where r)yn are the pseudorapidities of the jets. 

Since we will only consider unpolarized scattering, the correlation matrix element (Soft part) for the unpolarized 
hadron can be simply decomposed as [J, 0| 



fi{x,k 2 A 



M 



(9) 



where f\(x,k 2 A _) is the unpolarized TMD quark distribution, and hy (x, kj_) is the Boer-Mulders function describing 
the transverse spin distribution of a quark in an unpolarized hadron [2(]. 

In di-jet production, the transverse momenta of the particles in the di-jet are measured in order to probe the 
corresponding parton transverse momentum, and in this sense the azimuthal angle dependence of the di-jet is of 
most interest. The cross-section of unpolarized di-jet production contains both an azimuthal angle independent part 
and an azimuthal angle dependent part. For the azimuthal angle independent part, all possible partonic scattering 
subprocesses can contribute to the process, while only the unpolarized TMD distributions fi(x,Pj,) does. On the 
other hand, for the azimuthal angle dependent part, the Boer-Mulders function h± is involved, and since hi is a 
chiral-odd TMD function, only partonic processes that flip the helicity of the quark will contribute. The possible 
partonic processes are q(k%) +q(k 2 ) — ► q{Pi) +q(P2) and q(ki) +q(k 2 ) — > q(P%) +q{Pi) (or q'(Pi) + (f (P2)). Since the 
contribution of h^ 9 is power suppressed, we will not consider the gluon contribution to the azimuthal angle asymmetry 
in unpolarized di-jet production. Therefore we will specifically study the azimuthal angular dependence of di-jet in 
unpolarized Ha + Hb — > J q + J q + X and in Ha + Hb — > J q + Jq + X . However, since quark-channels make only a 
small contribution to jet rates (unless Pj_ is very large) this suppression might be compensated by the abundance of 
gluons and large scattering cross sections for gluons. 



A. Process H A + H B -> J q + J q + X 



For the process Ha + Hb — > J q + J q +X, the corresponding partonic subprocess is qq — > qq, and the lowest-order cut 
diagrams are shown in Fig.[TJ According to Eq. Q, one can therefore arrive at the angular independent cross-section 
given by: 
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FIG. 3: Partonic diagrams contributing to di-jet production Ha + Hb Jq + J q + X through gg — > qq. 
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da°[J q J q ] _ a 2 \N 2 -l (\ + {l-yf t 1 + y 2 \ (TV 2 - 1) 



d 2 P 1± d 2 P 2 ±dr ll dm ss \ 2N 2 \ y 2 (1 - y) 2 ) N 3 y{l - y) 

S 2 J d 2 k 1± d 2 k 2± S 2 (k 1± + k 2± - P 1± - P 2A _)fl {x u k 2 1± ) ft (x 2 ,k 2 2± ), (10) 



with s = P^_/(y(l — y)). The momentum fractions x\ and x 2 can be expressed as 



x x = ^(e" 1 + e" 2 ), X2 = ^( e -'n +e -^)_ (n) 
V s V s 



After integrating over P 2 ± and <pi, this leads to the expression 

d<T°[J g Jq] _ ir (N 2 -l (l + (l-y) 2 1 + y 2 



dPldr)xdr) 2 ss \ 2N 2 \ y 2 (1 - y) 2 

(A^ 2 - 1) 1 



N 3 y(l -y) 



}E e '/i(^)/i(^), (12) 



and the parton transverse momentum integration is de-convoluted by the 5-function. In the above equation we have 
used d 2 P 1± _ = \dPf x d<f>x « \dP\d^. 

Our purpose is to study the azimuthal angle dependence of unpolarized di-jet production. We recall that similar 
studies of angular dependence have been performed for the case of the unpolarized Drell-Yan process [III- Due 
to the correlation between the quark transverse spin and transverse momentum, the product of two Boer-Mulders 
functions x h j~ ca n give rise to large cos 2<f) angular asymmetries of dilepton production in unpolarized Drell-Yan 
processes [Til. l3l"l.l32T| . One can expect that the same mechanism can produce non-zero angular dependence of di-jets 
in unpolarized hadronic scattering. Actually a recent study (30j on photon-jet hadronic production shows that hi 
gives an angular dependence very similar to the cos 2ip asymmetry in Drell-Yan processes. We will follow this line to 
study the azimuthal angle dependence of unpolarized di-jet production. 

Since h x is a chiral-odd function which flips the helicity between the incoming and ongoing quarks, among the 
diagrams shown in Fig. [TJ only Fig. QJ: and Fig. [TJi can couple with h x and consequently contribute to the angular 
dependence of J q + J q production. Thus according to Eqs. ^ and ©, the contribution of hj^ to the cross-section of 
unpolarized J q + J q can be expressed as 



da A [J q J q ] 



{</(! - y)} E e « / d 2 k 1± d 2 k 2± S 2 (k 1± + k 2± - P x 



d 2 P 1± d 2 P 2± di ll df l2 M 2 sP[ 

k 1± ■ k 2± )(P 1± ■ P 2± ) - (k 1± ■ P 1± )(k 2± ■ P 2± ) - (k 1± ■ P 2± )(k 2± ■ P 1± ))hi q {x 1 ,k\ A _)hf q (x 2 ,fc 2 ± ),(13) 



where C qq is the color factor which will be determined below. This result is very similar to the one obtained in photon- 
jet production, and given in Ref. (30j . where a cos 20 angular dependence of photon-jet production is identified. Here 
we will use a weighting procedure on Eq. (|13[) . The advantage is that the integrations over k\± and k 2 ± can be 

deconvoluted. We find that using the weighting function W = -^-rp— - and integrating the cross-section in Eq. (fl3|) 
over P 2 ±, with 5(f> — n — ((f> 2 — (f>i), we can arrive at a simpler form: 

da™\J q J q ] f ^ id2p ^P 2 co S 5<f> 



dP 2 d m dr! 2 J "">"■"■ m 2 



da A [J q J q ] 



d 2 P 1± d 2 P 2± d m d V2 



S-L I 



(14) 



where h x (x) is the first fc^-moment of hi(x, k^) 



I 



K^(x)= I d'k^hf(x,ki) (15) 
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Eq. (fl4|) shows that the angular dependent part of the di-jet production cross-section can be obtained by a cos 6<p- 
moment of the total cross-section. 

Now we will calculate the color factor C qq appearing in Eq. (|14p . There are two sources which contribute to C qq . 

One is the color factor from the hard scattering subprocess. For the diagram in Fig. [TJ;, this color factor is — N am^ i 
which is the standard result of the generalized parton model. The other source is the special factor due to the 
presence of T-odd distributions, the so call gluonic pole factor 0, |2l[. The key ingredients corresponding to this 
factor are the initial- and final-state interactions (the gluon exchanges before and after the hard scattering), which 
lead to the gauge-links appearing in the definition of TMD T-odd distributions. The simplest example is that every 
T-odd distribution function will contribute a factor of +1 in SIDIS and —1 in Drell-Yan process, due to the final- 
and initial-state interactions, which correspond to the future- and past-pointing gauge links, respectively. In the case 
of hadron scattering the situation is more involved, since the initial- and final- state interactions are mixed, which 
makes the determination of gluonic factors not straightforward. Systematic study of the gauge-links appearin g in 
hadron scattering has been performed in Refs. [H, [H[. Here we will follow the approach presented in Refs. [l9l l2l| 
to calculate the gluonic pole factor coming from the product of two . The general quark correlator with gauge link 
U is given by: 



(2tt)2- 



This gauge-link is process dependent and also depends on the hard subprocess. As we mentioned before, the subprocess 
shown in Fig. [1}; can couple with and contributes to the angular dependence of di-jet in hadron scattering. The 
gauge link for the quarks which participate in this hard partonic process (the interference of t and u channel) can be 
written as [l9j : 

u m ~ u m ~ N i_ x N 01 

_^±1 W RI W H (17) 

where 

W P= W Ht W [+] (18) 

is the gauge link loop, and and U\~\ are the path-ordered future and past pointing Wilson lines. The two terms 
appearing in Eq. (|17p correspond to the contributions from different color flows of the qq — ► qq process in the tu* 
channel. The convolution of the soft parts in the decomposition of two color flows is therefore 

$M <g) $I wt ] 2N2 $[(■=■)+! (g> 

iV 2 — 1 



N z + 1 
~N 2 - 1 



$[□+] g, (19) 



where (□)+ and □+ are the shortcuts of the basic gauge links Tr ^. — -l/W and i/W, respectively. 

Here we will consider only the T-odd part of the contribution, the so call gluo nic-pole contribution, since we are 
interested in the contribution of products of two h ± . As shown in Ref. [l^, |21| , after weighting by the transverse 
momentum, the T-odd part of each basic soft correlator will appear as the T-odd distribution function (in our case 

the T-odd functions is ft-i" ) multiplied by a prefactor. For example, the correlator <£>[( n )+l will give a factor 1, while 
will contribute a factor 3. Therefore, the factor contributed by the product of two is: 

2N 2 N 2 + 1 

x 1 x 1 x 3 x 3 



N 2 -\ N 2 - 1 

7A^ 2 + 9 
N 2 -l 



(20) 



We emphasize that this factor is the result coming from the T-odd distributions due to their complicated gauge-link 
structure, which cannot be predicted by the generalized parton model. After considering the color factor from the 
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hard process, we get C qq as: 



_ TV — 1 / IN 2 + 9 
99 ~ 47V 3 " X \~ TV 2 - 1 
7iV 2 + 9 



(21) 



47V 3 

Therefore the final result of the cos ^0-moment shown in Eq. (Tj"4"| is 

rfg w [J g J q ] A-wal IN 2 + 9 
dP 2 d Vl d m ~ sP 2 N 3 y ^~ y > 

x XX (1) '(zi)^ (1), (z 2 ). (22) 
We then can define a subprocess cos asymmetry in if 4 + Hb J q + J q + X processes as 



^"[JgJj/dPjrf^ _ 4(77V 2 + 9) J2 q hj (1)q (x 1 )ht (1)q (x 2 ) 

The main result given by this asymmetry is that there is an additional factor — ^ jjp , as shown in (f20|) . apart from 
the one obtained in the standard generalized parton model. The size of this gluonic pole factor is -9 for N — 3, which 
is quite large and will greatly enhance the asymmetry. Moreover, comparing to the result from generalized parton 
model, the sign of the asymmetry is reversed by this gluonic pole factor. 



B. Process Ha + Hb — > Jq + J- q + X 



Now we will consider the process Ha + Hb — > J q + J q + X , where the partonic subprocesses contributing include 
qq — > qq and g q — > q q, as shown in Figs. [2]and[3l respectively. The diagram in [2^, is also present for qq — > q' q 1 , which 
will contribute to J q + J q di-jet production. As before, the angular independent cross-section of Ha+Hb — * J q + J q +X 
process can be readily written as 



*oy. H _ 4 t d , kl±d > kl±5Hkl± + k2± _ Pl± _ p 2± j ; f l±£-^ + (1 _ yf + / 



d 2 P 1± d 2 P 2 ±d m dr]2 ss J ^ ^ v ^ ^ \\ N 2 V 2 2/ 

, (iV 2 -l)(l-y) 



iV 3 y 



J2(f! (xi,k 2 1± ) ft (x 2 ,k 2 2± ) + ( Xl ~x 2 )) 



1 " 2 ' " " ]2 V y 2 (i-y) 2 )f 9 i(xuk 2 1± )fl(x 2 ,k* ± )\, (24) 



,27V y(l -y) iV 2 - 1 
and the corresponding P2J_-integrated cross-section is given by 



+ — JLw'WiW 

- *»)) + (^^r^p - iv^ivd - y)) A 9 (*0/f (*»)} , (25) 

Again, since /i^ is a chiral-odd function which flips the helicity of the initial and final partons, the hard diagrams 
that can contribute to the angular asymmetry in Ha+Hb — > J q + J q are given in Fig. [2^, Fig. [2]: and Fig.[2ji. Following 
a similar procedure as that used in the previous subsection, we can give the cos <5</>-moment of Ha + Hb — > J q + J q 



da° [J q Jq] na 2 
dP^di]idr/2 ss 



N 2 -l / l + (l-y) 2 
TV 2 v. 2y 2 



1 - 2y + 2y 
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FIG. 4: The subprocess cos 8<f> asymmetry R qq in di-jet production Ha + Hb — > J q + J q + X at ^J~s — 200GeV, for Pj_ = 10 GeV 
(left column) and 20 GeV (right column). The dashed line and solid lines are the results from taking Boer-Mulders functions 
of Set I and II, respectively. The upper and lower panels show the asymmetry vs y and 771 , respectively. 
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FIG. 5: The P^-dependent subprocess cos Sep asymmetry R qq in di-jet production Ha + Hb — » J q + J q + X at y/s = 200GeV 
calculated from the Boer- Mulders functions in Set II. 



using the weighting function W = P± ^ S ^ , as follows 



dP±dr)idr)2 J d 2 Pi±d 2 P 2 ±dr)idr]2 

' i,;Tn? (C qq y*(l-yf + C > qq y\l-y)) 



x 



sPl 

(26) 

The factors C qq and C' qq can be calculated by combining the contribution of the hard subprocess and the gluonic pole 
factor. The term with C qq is the result coming from the hard process in Fig. [2^,, and the color factor of this process is 

■ The term with C' qq is the result from the hard processes Fig. [2k and Fig.[3H. The color factor of these processes 

is — ^jyis 1 ■ The gluonic pole factors for those processes can be calculated from the corresponding gauge-links. For the 
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diagrams in Fig. [2^, c and d, the gauge-links are equivalent, and given by: 

u [ss'] u [st>] = u [ts>] N 2 [(D)+] 

t^qq ^qq ^qq jy 2 _ ^ 



N 2 - 1 

The convolution of the soft parts in the decomposition of two color flows is therefore 

„ 2N2 $[( d )+] $[( D )+ t ] 
A'' 2 — 1 



1 $H. (27) 



' $Hg, (28) 



iV 2 - 1 

Therefore the factor contributed by the product of hj^ and is calculated as: 

iV 2 1 

— = xlxl = x (-1) x (-1) = 1. (29) 

N 2 - 1 iV 2 - 1 V ' V ' V ' 

Thus we conclude that C qq — N 4 ^i and C' q = — ^^g 1 - Unlike in the case of J 9 + J 9 production, this result 
shows that in the case of J q + J q production there is no color factor enhancement for the azimuthal angle dependent 
cross-section, being therefore equivalent to the generalized parton model result. Thus we can give the subprocess 
cos 8(f) asymmetry of J q + J q as the ratio of Eq. (|2l)j) and ([2!>]) : 

_ da w [J q J q }/dP 2 dr h d m = 

Mv) E q ( /v (1) ' 9 (^i) hf'^( X2 ) + ( Xl «-> x 2 )) 

(3D) 



B{y) Eq(f? (si) fl (* 2 ) + (xi - 12)) + C(y)ff ( Xl )ff (x 2 ) 



where 



A(y) = 4 ( y (l -*)-■£), (31) 



1 + (1 - yf 

1 (1-y) 2 



B(y) = - ' v 2 2 " y +(i-y) 2 + y 2 



iV y 

TV 2 - 1 y 2 + (1 - y) 2 1 



(32) 



c ^ - ' 2^3 %(i-,r -N y{i - y) ^ (33) 



III. NUMERICAL RESULTS 



In the previous section we have calculated the theoretical result for the azimuthal angle dependence of di-jet 
production in unpolarized hadron scattering, and we showed that the azimuthal angle dependent part of the cross- 
section can be separated from the azimuthal angle independent part by taking a cos (50- moment. We defined the 
cos 5(f) asymmetry as the ratio of cos 50-moment and the azimuthal angle independent cross-section. It is interesting 
to study wether the asymmetry could be accessed in hadron colliders. The Relativistic Heavy Ion Collider (RHIC) at 
BNL is running the polarized di-jet production process p^ +p — ► J\ + J 2 + X [33[ , with the main purpose of measuring 
the transverse single spin asymmetry. The process can also be used to access the unpolarized di-jet production by 
averaging the spin-dependent data. We therefore calculate the subprocess asymmetry R qq given in Eq. (f2"3"|) at RHIC 
energy [yfs = 200 GeV), and show the result in Fig. |U In the calculation we take N = 3 and apply two sets of 
Boer-Mulders functions. Set I is the Boer-Mulders functions which saturate the positivity bound: 

hf^(x) < (34) 

with (|fej_|) = 0.44 GeV 34]. Therefore this Set Boer-Mulders function will give the maximum bound of the asymme- 
tries. In Set II we apply the Boer-Mulders functions which are extracted in Ref. [351 ] from unpolarized p + D Drell-Yan 
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FIG. 6: The subprocess cos Sep asymmetry R q q in di-jet production Ha + Hb — > Jq + Jq + X at = 200 GeV for P± = 20 
GeV. The dashed line and solid lines are the results from Boer-Mulders functions of Set I and II, respectively. 



data 0] measured by the FNAL-E866/NuSea Collaboration. As shown in that paper, the Boer-Mulders functions 
in Set II can successfully reproduce the cos 2c/> asymmetry in the unpolarized p + D Drell-Yan process. So far the 
information on these functions is limited, therefore we adopt these two sets of Boer-Mulders functions to perform 
estimates on the cos 8<j) asymmetries. For the unpolarized parton distributions appearing in the denominators of 
Eqs. (23]) and ([30]) adopt the MRST2001 (LO set) parametrization The results for R qq from sets I and II are 

plotted with dashed and solid lines, respectively. The upper panel of Fig. 2] shows the subprocess asymmetries vs y 
at P± = 10 GeV (left) and P± = 20 GeV (right) respectively, while the lower panel shows the subprocess asymmetries 
vs rji at same P±. In Fig. [5] we also show R qq vs P± which is calculated from the Boer-Mulders functions in Set II. 
The result from Set I, which can be understood as the upper limit of the asymmetry, can be more than 10% in size. 
The result from Set II is about several percent, which is still quite sizable. The main reason of this result is the 
enhancement factor of -9 on the asymmetry from T-odd distributions in hadron scattering. 

We also calculate the subprocess cos 5(f) asymmetry of J q + J q production given in Eq. (j3"0)) at \fs — 200 GeV, and 
show y-dependent curves at P± = 20 GeV and r\\ = 0.2 in Fig. [S] Our calculation indicates a subprocess smaller 
asymmetry for J q + J q production than that for J q + J q production. This can be understood from the fact that there 
is no color factor enhancement for the azimuthal angle dependent cross-section of J q + J q production, and there is a 
large gluon-gluon scattering contribution in the denominator of Eq. (|30|) . 

In a practical measurement of di-jet production, the flavor of the jets usually is not measured. Hence the subprocess 
asymmetries shown in Figs. HI E] and [5] are not measurable. The purpose to present those figures is to show how much 
the contribution coming from each subprocess is. Therefore we calculate the cos 5<fi asymmetry of any di-jet production 
Rjj, which is measurable and defined as 

R _ Ea, b=q , q , g da W [JaJ b }/dPld m d m 

" Ea, b = q , q , g d<y°[JaJ b }/dPld m dr, 2 - 1 j 

In this paper we assume that the partonic processes that contribute to the azimuthal angular dependence of di-jet 
production are qq — > qq and qq — > qq. Therefore the numerator in Eq. (|35[) is the sum of Eqs. (|22p and (|26p . while the 
denominator contains the contributions from Eqs. (fl2|) . ([25]) and from all other possible processes, such as qq' — > qq', 
qg — > qg, and so on. In Fig. [7] we show by thin lines the asymmetry Rjj at RHIC for P± — 20 GeV vs y, from the 
two sets of Boer- Mulders functions (solid and dashed lines for Set I and Set II respectively). As a comparison, we 
also show by thick lines the result of Rjj from the conventional generalized parton model (GPM) calculation where 
gluonic pole contributions are not considered. In Fig. [7]d, we show the asymmetry Rjj at RHIC for Pi = 10 GeV 
(thick lines) and P± = 20 GeV (thin lines) vs r)i + r]2- In the calculations we restrict — 1 < r\xn < 2 to consider the 
pseudorapidity acceptance of jet at RHIC. The curves given in Fig. [7K indicate a sign reversal of the asymmetries 
Rjj between the results from the two different approaches. This is due to the additional factor —9 in the azimuthal 
angular dependent part of Ha + Hb — > J q + J q + X , coming from the multiple initial- and final-state interactions. We 
remind that a similar sign reversal also happens in the SSA in hadronic photon-jet production (27j . which is due to 
a similar effect. Therefore the measurement of the cos 8<f> asymmetry Rjj in unpolarized di-jet production at RHIC, 
especially the enhancement of the size and the sign reverse of the asymmetry, will provide further tests on the role 
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FIG. 7: (a), The y-dependent cos 8<j> asymmetry Rjj in di-jet production Ha + Hb — » Ji + J2 + X at RHIC for P± — 20 GeV. 
The dashed line and solid lines are the results from Boer-Mulders functions of Set I and II, respectively. We also show by thick 
lines the result of Rjj from the conventional generalized parton model (GPM) calculation where gluonic pole contributions are 
not considered, (b), The cos 5<j> asymmetry Rjj vs 771 + 772 in di-jet production Ha + Hb — > Ji + J% + X at RHIC for P± = 10 
GeV (thick lines) and P± — 20 GeV (thin lines) from Boer-Mulders functions of Set I (thin lines) and II (thick lines). 



of initial- and final-state interactions as QCD dynamics of T-odd distributions in hadron scattering. We would like 
point out that in calculating Rjj, we only considered the lowest order contributions to parton scattering. Other 
contributions may affect the size of the asymmetry, such as NLO corrections for the denominator of Eq. (|35|) . At large 
P±_ they are further enhanced by threshold logarithms which are more pronounced for gluonic channels. Therefore 
our results are estimates and Rjj can be smaller. However we expect that those contributions will not change the 
sign of the asymmetry. 

IV. SUMMARY 

In summary, we performed a study of the azimuthal angular dependence of back-to-back di-jet production in 
unpolarizcd hadron scattering, focusing our attention on the process H a + Hb — * J q + J q + X, where two quark jets 
have the same flavor, and on the process H a + Hb — > Jq + Jq + X, where the di-jet comes from a quark-antiquark 
pair. We find that in both processes there are an azimuthal angle dependent cross-section generated by the product of 
two hi from each incident hadron. Using a weighting function cos 6(f>, we can separate the azimuthal angle dependent 
cross-section from the azimuthal angle independent one. We defined the cos 8<j> asymmetry of di-jet production as the 
ratio between the cos (50-moment and the azimuthal angle independent cross-section. In the case of J q + J q production, 
due to the multiple initial- and final-state interactions within hadron scattering, there is a color factor enhancement 
(which is -9 for N = 3) in the asymmetry compared to the result from the standard generalized parton model. No 
such enhancement appears in the case of J q + J q production. We estimate the cos 8<f) asymmetry of di-jet production 
at the RHIC energy, and find that the subprocess asymmetry of J q + J q production can reach up to 10% in maximum 
due to the large color factor enhancement. In contrast, the subprocess asymmetry of J q + J q production is smaller. 
We further find that the color factor enhancement in J q + J q production can increase the size and reverse the sign of 
the total cos 8<fi asymmetry Rjj in di-jet production. Therefore it is feasible to perform the measurement on the di-jet 
production in unpolarized hadron scattering at RHIC, which can identify the color factor enhancement on the cos<50 
asymmetry in J q + J q production. Furthermore, the study will provide the opportunity for a better understanding on 
the role of initial- and final-state interactions in hadron scattering and other processes. 
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